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Abstract 



It is well known that the Pauli principle plays a substantial role 
at low energies because the phonon operators are not ideal boson op- 

■ erators. Calculating the exact commutators between the quasiparticle 
£C) \ and phonon operators one can take into account the Pauli principle 

■ corrections. Besides the ground state correlations due to the quasi- 
particle interaction in the ground state influence the single particle 
fragmentation as well. In this paper, we generalize the basic QPM 

\q \ equations to account for both mentioned effects. As an illustration of 

■ our approach, calculations on the structure of the low-lying states in 

131 Ba have been performed. 

1 Introduction 

In the forthcoming period there will be an increasing activity in the domain 
of unstable nuclei studies due to the start of operation at several major 
facilities. A theoretical investigation of odd nuclei far from stability demands 
to include the ground state correlation effects. 

The quasiparticle-phonon nuclear model (QPM) is widely used for 
the description of the energies and fragmentation of nuclear excitations. 
The different versions of the QPM equations for odd spherical nuclei are 
given in [21 |3J |1] . It has been shown j^l [3] that the Pauli principle plays 
a substantial role at low energies, but the ground state correlation effects 
were not taken into account in these calculations. From the other side, the 
ground state correlations influence the single particle fragmentation |Sj as 
they shift the strength to higher excitation energies. 
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In this paper, we generalize the basic QPM equations to account for 
both mentioned effects. We treat long-range ground state correlations by 
including backward-going quasiparticle-phonon vertices using the equation 
of motion method 7 with explicitly accounting for the Pauli principle. As 
an illustration of our approach, calculations of the structure of the low-lying 
states in 131 Ba have been performed. 

2 Basic formulae 

We employ the QPM-Hamiltonian including an average nuclear field de- 
scribed by the Woods-Saxon potential, pairing interactions, isoscalar particle- 
hole residual forces in separable form with the Bohr-Mottelson radial depen- 
dence 8_: 



(iiP) , , 

H = EOi " X r)a) m a ]m - -G^ : {P%Ptf : -- £ «W : (M^,) :}. (1) 

t jm Xfi 

The single-particle states are specified by the quantum numbers (Jm); 
Ej are the single-particle energies; A T is the chemical potential; and 
k( a ) are the strengths in the p-p and in the p-h channel, respectively. The 
sum goes over protons(p) and neutrons(n) independently and the notation 
r = {n,p} is used. The pair creation and the multipole operators entering 
the normal products in Q are defined as follows: 

jm 

M i = ^X+T ^ fjjhjmj'm | AM)ot m a J ./ m / 1 

jj'mm' 

where /y are the single particle radial matrix elements of the residual 
forces. 

In what follows we work in quasiparticle representation, defined by the 
canonical Bogoliubov transformation: 



+ , / -I \ j—m 
U j a Jm + (- 1 ) VjOtj 



"jm - "J^jm ^ V "j«j-m- 

The Hamiltonian can be represented in terms of bifermion quasiparticle 
operators (and their conjugate ones): 
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The phonon creation operators are denned in the two-quasiparticle space 
in a standard fashion: 

Qt* = \ Eti$ A+ ^'-> A ^ ~ (- 1 ) A ~ / V& - m)}, 

. ./ 

33 

where the index A = 0, 1, 2, 3, ... denotes multipolarity and \x is its z-projection 
in the laboratory system. The normalization of the one-phonon states reads: 

In terms of quasiparticles and phonons the Hamiltonian is rewritten 
H = ho + h pp + Iiqq + tiQB, 

ho + hpp = £ j a ^m a 3 m ' 

jm 

h QQ = -\ E ^(^') + (-) A "" oa-^) + h a+|,< V)> 

A/m 

^ = "27! S ^r( J i'A i )((-) A ^Q+ Mi + Q A ^)i?(i/;A-^ + / l . C ., 



where 



„(-) f w 

33 33 



* A * = E 



2 , ,2 ' 



" J (1/ - ^0 



3 



with v-J = UjUji — VjVji and u--, = UfVj + UjVj>. Here and further we use 
the notation 7r,- = w (2j + 1). 

The model wave function of an odd spherical nucleus is taken in the form 

* u (JM) = Oj Mu \), (2) 

where 

+ jmu = Cj»<* + JM+Y. D T (Jv)P^{JM)-E Ju aj M -Y, F j\ J v) P MJM\ 

i i 

(3) 

with 

pU jm ) = 

and ~ stands for time conjugate according to the convention: Pj\i(JM) = 
(-1) J - M PMJ-M). 

We apply the equation of motion method to the excitation operator ©: 

(\{60j M „,H,0+ Mu }\)=r)j v (\{60jM,0+ M }\). 

Following the linearization procedure [7j, at the final state of calculation 
of the matrix elements, we consider the ground state to be a vacuum state 
for both operators oljm and Qx^i- 

In all calculations the exact commutation relations between the quasi- 
particle and phonon operators are considered: 

[(XjmiQjMv) = Yl (W m ' I JM )^jp a f m'- 

j'm' 

The normalization condition of the wave function reads 
(\{Oj M u,Oj M J\) = C% + E% + Y}Df(Jv)f + Y}F}\Jv)?+ 

j\i jXi 

+ Yl [D^(Ju)Df(Ju) + F^(Ju)Fp' i \ju)}^j(j\i\j'X'i') = l. 

jXij'X'i' 

The equation of motion leads to the following system of linear equations 
for each state with quantum numbers JM : 



a jm Q 



X/ii 



.JM 
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/ ej V {Jj'X'i') G{J) -W(Jj'X'i') \ 

V(JjXi) Kj(jXi\j'Xi') -W(JjXi) Sj(jXi\j'X'i') 



G{J) 



-W {Jj'X'i' 



-ej 



-V (Jj'X'i') 



\-W(JjXi) Sj(j'X'i'\jXi) -V(JjXi) -Kj(jXi\j'Xi') J 



D 3 l ( J ") + E Df'{Jy)^j{jXi\j^i') 



(C Jv 

Df'{Ju) 

—Eju 
\-Fp'{Ju)J 



j'X'i 1 



-Ej u 



-Fftju) - Fp'(Jv)J?j(jXi\j'XY 

j'X'i' 

The average value of H over the wave functions (J2J) is 



V 



(4) 



< Vt(JM)HW v (JM) >= 

= [Cj u -E%]ej+ £ (Df(Ju)Df\ju) + F^{Ju)Ff{Jy))Ij(jXi\j'X'i') 

jXij'X'i' 

+ 2j2(Cj u D^(Ju) - E Ju F 3 Xi (Ju))V(JjXi) + 2Y,{CjuFf{Jv) + E Ju Df{Ju))W{JjXi) 

jXi jXi 

-2C Jv Ej u Gj-2 D^Fp'Sj(jXi\j'X'i'). 

jXij'X'i' 

We give the explicit expressions for the quantities entering the above 
formulas with short comments 



J?j(jXi\j'X'i') = vtavtv E^iViu' 

h 



f h A 

j J X' ' 



Gj = (\{S5 M , [H,a+ M ]}\) = V2^^r(JjA,)^, 



Xij 



V(JjXi) = (\{[a JM ,H],P+ M (JM)} 



V2 



r ( JjM) 



-4 E U^i'M) (jAi|/AV)) r (Jfx'i') 

V 2 j'X'i' 
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As a result of the application of the equation of motion method, the 
matrix elements V (JjXi) between quasiparticle and quasiparticle(^)phonon 
states differ by the ones obtained earlier [2j by an additive containing (jXi\j'X'i') 



W(JjXi) = (\{[a+ M ,H],P+ Xi (JM)} 



I TO 



£ ^ (X'i'i") — ^:^j{jXi\j'X'i'')-^f^j[jXi\j'>!i') 



\'j'i'i"T 



TTJ L 



The matrix elements W (JjXi) appear after the inclusion of the backward- 
going terms in the operator @ and they present a central issue of this work. 

Sj {jXi\j'X'i') = (\{P+ Xi (JM) , [H,P+ yi , (JM)]]\) = -Gf&j , 



I j (jXilfX'i') = (|{p jAi (JM) , [H,P+ xli , (JM)] }|) 



Ij(jXi\j > X'i > )+I J (j > X'i'\jXi) = 

= 25 jj >5xx'S iil (uj X i + £j) + ££j (j'X'i'\jXi) [syj + + ^Xi) - {jXi\fX'i') , 
Kj(jXi\j'Xi') = X - [Ij (jXi\j'X'i')+Ij (j'X'i'\jXi)] 

(jXi\j'X'i') = i £ K> (A*ii) (f X'i'\jXh) + < (X'hi') ££j {jXi\j'X'h)] + 
+ 7 E < ( Al "') (J Xi \nXik) %j (/AYjiiAaia) + ifj (/AVtfiAaii) if; tfAi|jiA a 



Aiiii2j'iTo 



The quantities .i?j(jAi| j'AV) , 5j(jAi| j'X'i') and (jXi\j' X'i') vanish 
if the Pauli principle is not respected. 
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3 Approximations 



3.1 General 

As has been shown in [2] «Sfj are alternating quantities and their diagonal 
values are much greater than the nondiagonal ones. This is natural from 
the physical point of view as the Pauli principle is violated most probably 
in the configurations formed by identical quasiparticles. The same applies 
for the new quantities SFj. 

J?j(jXi\j'X'i') = SfiJjXWwSwSw, 

^j(j\i\j'x'i') = 3r{jjxi)s jf s xxl d w , 

where 

In this approximation 

v (JjXi) = ~{i + sf(jj\i) + sr{JjXi))T (jjXi) . 

The vertice V is renormalized by the factor (1+jSf (JjAi) + 3? (JjXi)). In 
configurations with strong Pauli principle violation the quantities ££(JjXi) 
go to —1 and the quantities S?(JjXi) go to 0. The role of the term 2? (JjXi) 
in the renormalization becomes more important as the phonon collectiveness 
increases. 

W (JjXi) = + J?(JjXi) - J?(jJXi)) £ < (An') 

4 7TJ z — ' J 

i'ro 

As with the vertices V, the vertices W are renormalized now by the 
factor (1 + J£(JjXi) — ££(jJX%)) and again in configurations with strong 
Pauli principle violation the quantities ^£(JjXi) go to —1 and the quantities 
Jzf(jJXi) go to 0. The results for V and W show that configurations with 
J? (JjXi) close to —1 must be excluded from the configuration space. 

Kj (jXij'X'i') = S jf S XX 'Sii' (1 + ^(JjXi)) (u M + ej-M (JjXi)) , 



(j 3 xi)-^ x ^ Xi) 



l + 3f(JjXi)' 



The quantities Si (JjXi) play a very important role as they shift the 
values of the poles and this shift depends on the extent of the Pauli principle 
violation [2]. 

Neglecting Gj and Sj(jXi\j'X'i'), we arrive at the system of equations 











+ 



Mu Mia 
M 2 i M 22 



-Ej u 



-Ej u 



where 



V 2 (JjM) 



+ 



W 2 (JjM) 



j£ (1 + ££ (JjXi)) \r] Jv - {com + Sj-M {JjM)) r} Jv + uj X i + Ej-& (JjM) J ' 

(5) 



M 22 = £- 



jXi 



W 2 (JjXi) 



+ 



V 2 (JjXi) 



(1 + 3? (JjXi)) V??j^ - (ujm + ej-M (JjXi)) 7] Jv + u Xi + Ej-M (JjXi) J ' 



Mi 2 = M 2i = 

V(JjXi)W(JjXi) ( 1 



a 



^ (l+jSf(JjAi)) \r} Ju + uJxi + £j -&(JjXi) r} Jv - (u Xi + £j - St (JjXi)) /' 
leading to the equation 

M 12 M 21 = (ej + Mn - t?^) (M 22 - ej - n Jv ). (6) 

3.2 Limit cases and analisys 

The equation © can be approximated by the following one 

M 2 

ej -^"- M "- fe-(%-M ll)r (7> 

Therefore, neglecting the backward amplitudes, i.e. setting W (JjXi) = 
0, (0) immediately reduces to the secular equation obtained earlier [2] : 

F _ n = v ^ 2 

J ^ |- (( WA< + £j - & (JjXi)) - VJu ) (1 + J2P ( JjAi)) ' 1 j 

The significant difference for the solution rjj of the equation (J7J) as com- 
pared to the equation comes from the second term in the r.h.s. of the 
expression (jSJ which contributes to a shift of the first solution of equation 
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(|7|) to higher energies. The second term of the r.h.s. of equation also 
contributes in the same direction but to a much smaller extent. The shift 
in energy becomes larger as the interaction between the quasiparticles and 
phonons increases. A critical value for the interaction exists as in (jHJ) but 
now due to the second type of terms in M%% an increase in the interaction 
leads to a shift of the first solution towards the pole as contrary to the case 
neglecting backward amplitudes where the solution moves in the opposite 
direction. 

The inclusion of the matrix element Gj in the system of equations can 
be accounted as a change of M\i to M12 + Gj but the resulting shift in the 
first solution turns out to be negligible. 

4 Numerical results 

In order to give a qualitative picture of the effects on the structure of the 
low-lying states imposed by the backward-going amplitudes, numerical cal- 
culations for 131 Ba were performed. This isotope belongs to the transitional 
region where the anharmonic effects play a gradually increasing role at low 
and mainly at intermediate energies and therefore the results presented in 
this section may lack some accuracy because the wave function (J2J) does not 
contain configurations to account for these effects. The parameters of the 
Woods-Saxon potential are as follows: 



A 


N,Z 


ro, fm 


V , MeV 


k, fm 2 


a, fm 1 


G N>Z , MeV 


127 


N = 74 


1.280 


43.40 


0.413 


1.613 


0.124 


127 


Z = 53 


1.240 


59.72 


0.350 


1.587 


0.130 



Table 1. Parameters of the Woods-Saxon potential for 131 Ba 



Our study shows that in realistic calculations one must include phonons 
with A = 2, 3, 4, 5. The strength parameters are adjusted so that the 
odd energy spectrum of the low-lying states is reasonably close to the ex- 
perimental values. As a result, the energy of the first quadrupole state of 
130 Ba has a value that is much higher than the experimental one. 

Solving the systems of equations (jljl. one can find the structure of the 
wave functions (J2J) and the energies of the excited states. Working in a di- 
agonal approximation for S£j and 3Tj this system reduces to a generalized 
eigenvalue problem. In fig.l a comparison between the experimental values 
of the energies and the theoretical calculations is presented. We restrict the 
calculation to the six neutron states 1/2+, 3/2 + , 9/2~, 11/2", 5/2 + , 7/2 + . 
The level ordering presented in the third column on this figure generally 
agrees with the one obtained in [3]. The results clearly support the conclu- 
sion following from as the first solutions obtained after the inclusion of 
the backward-going terms become closer to the first poles and consequently 
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1/2+ 






2115/2 



Bxperiment 



Quasiparticles 



7/2+ 


2097 


9/2- 




3/2+ 


1635 


5/2+ 


13(19 




/HsT 


3/2+ 


141 1 


7/2+ 


993 


5/2+ 


744 


11/2- 


150 


9/T\ 






V if 


Without Backward Amplitudes 



3/2+ 


1071 


9/2- 


996 




839 






3/2-; 


g 786 






7/2+ 


505 


11/2- 


287 


9/2- 


187 


3/2+ 


108 


1/2+ 






With Backward Amplitudes 



Low-lying energy spectrum of 131 Ba (in KeV). The first column gives the experimental values 0, 
the second is the unperturbed lqp spectrum, the third gives the levels resulting from the solution 
of JSJ, the fourth is the full calculation with backward amplitudes and Pauli principle corrections. 



closer to the second solution thus significantly reducing the gap between the 
first l/2 + and the second l/2 + states as well as between the first 3/2 + and 
the second 3/2 + states. The intruder state 9/2~ deserves a special attention. 
The wave function of this state is practically a pure quasiparticle(^)phonon 
state with a structure [l/in/2 2^]g/2- • The significant reduction of the 
energy of this state is due to the Pauli principle correction, the inclusion of 
which is essential for the correct ordering of the first several levels. Along 
with the experimental energies our calculations provide a reasonable de- 
scription of the spectroscopic factors for the (d,p)-reactions: for the states 
1/2+ and 3/2+ having experimentally measured values of the spectroscopic 
factors of 0.53 and 1.03 respectively, our calculations give 0.46 and 1.36 . 



5 Conclusion 

The comparison between theoretical calculations and experimental data for 
131 Ba has shown that in order to describe the structure of the low- lying 
states in odd-mass nuclei far from the magic numbers one needs to take 
into account the Pauli principle and the ground state correlations effects 
simultaneously. Calculations for other Ba isotopes are in progress now. To 
improve this approach a self-consistent description of the mean field with 
more realistic effective nucleon-nucleon forces is desirable. 
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